A model for ultrarelativistic spherically symmetric 
Pre- Hawking radiating gravitational collapse 



Miguel Pinol Ribas 
March 3, 2013 

Abstract 

In this paper we present a simple but non-simplistic model of gravitational col- 
lapse with thermal emission of pre-Hawking radiation. We apply Einstein equations 
to a time-dependant spherically symmetric metric and an ultrarelativistic stress- 
energy tensor. In our model, particles either radially approach the center of the star 
as collapsing matter, or radially flee from it. 

1 The Model 

The more general metric for a spherically symmetric collapse presents the following 
structure: 

dr 2 = e v dt 2 - e x dr 2 - r 2 dtt 2 , (1) 

where v and A ara functions of radius r and time t. It will be useful to define also a 
function 0, so that 

u = -X-2(f) (2) 

and 

dr 2 =e- X - 2 Ut 2 -e X dr 2 -r 2 dn 2 , (3) 

as it may be demonstrated [Lan] that <f> > all over to space (it is strictly equal to 
zero in the "outer space", outside the edge of the collapsing star, where the metric is 
identical to Schwarzschild's one). In addition, (j)' < all over to space (we are denoting 
with x 1 derivation in respect to r, as we will denote with x derivation in respect to t all 
over the paper). 

1.1 The Einstein equations 

The Einstein equations for a metric of this kind are these ones: 
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SttTI = -e- x (- + \) +1 (4) 
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SttTq 1 = -e" A ^. (7) 
The other components vanish identically. 



1.2 The stress-momentum tensor 

The stress-momentum tensor for a perfect fluid takes a simple form, 

T$= g/35 (p + p)u a u s -ri%p, (8) 

where p corresponds to the total energy density and p to the pressure. For an 
ultrarelativistic fluid, p is proportional to p: 



p = up (9) 

(concretely, it may be demonstrated that uj = |, but we will keep on writing the constant 
of proportionality along the paper). 



1.3 The four-velocity in our model 

In our model, all particles (or " elements of mass" , if you prefer) describe a purely radial 
trajectory (inwards, or outwards), so that the only non-null components of four- velocity 
are u° (which we are going to call 7, in order to obtain more beautiful expressions) and 



u 



From the metric, we can deduce the relationship between both components: 



e A «) 2 , 



(10) 
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(where the subindex p stands for "particles"). 

For infalling particles, u r p = — Up ; for outgoing particles, Up = + 
if we call \ the proportion of infalling particles, so that 1— d be the proportion of outgoing 
particles, the mean "r- velocity" u 1 will be the following one: 



. Consequently, 
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1.4 The stress-energy tensor components 

The stress-energy tensor components are, with all the considerations that we have taken 
until now, the following ones: 

T ° = e~ A - 2 V (p + p) - p = e - A ~ 2 <V (1 + u,) p - up, (13) 



Tl = -e x ~/u 1 (p + p)-p=(l- 2x)e- < V V 1 ~ e A + 2< ?V 2 (1 + u) p - up, (14) 

T* = Tl = ~P = -up, (15) 

= e-^V 1 (P + P) = "(1 " 2x)e- 2A - 3 V\/l - e A + 2 *7~ 2 (1 + u) p. (16) 
1.5 The ultrarelativistic limit 

In the ultrarelativistic limit, 7 >> 1, so that we may ignore the terms with lesser powers 
of 7 in front of those in the higher ones: 

T «e- A - 2 V(l + ^)/>, (17) 
Tl^(l-2 X )e-^ 2 (l + u)p, (18) 

T 1 «-(l-2 X )e- 2A - 3 V(l + ^)p. (19) 

Therefore: 

/ 9«e A+2 V 2 (l+^)^ 1 T °, (20) 
Tl « (1 - 2 X )e A +*T °, (21) 
T 1 «-(l-2 X ) e - A ^T °. (22) 

2 Solution of equations 

2.1 Some considerations on the phases of collapse, in the light of pre-Hawking radi- 
ation 

In 2006, Vachaspati, Stojkovic and Krauss demonstrated that collapsing stars emit pre- 
Hawking radiation [Vac]. While there exist some difference between both, pre-Hawking 
radiation spectrum results to be roughly proportional of Hawking radiation of black 
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holes. Thus, we are going to use a Hawking radiation-like expression for pre-Hawking 
radition in our model of gravitational col-lapse: 



—k 

m p -H ~ —2, (23) 

where A; is a constant of proportionality, and m denotes the inner "mass" which is 
lost as "radiation" . In our model, we are not only going to consider the "global" emission 
of the collapsing star but also that of the "inner layers" towards the outer ones (and we 
will assume the same expression for that). 

If we have into account the relationship between the function A and the "mass", 

e A = 1 " — , (24) 
r 

we may straightforwardly deduce that 

e-*A = (25) 
r 

A = (26) 
r 

On the other hand, from equations (7) and (22) we may obtain another expression 
for A: 

A = -re A 8vrT 1 = r(l - 2x)e^87rT °, (27) 

which we may split into two terms in order to make explicit the contribution of 
infalling and outgoing fluxes: 

A = r(l — x)e~^7rT - r X e^87rT °. (28) 

Consequently, we may logically identify the second term in the previous equation 
with the pre-Hawking variation of A which we may obtain from equations (23) and (26): 

■ 2e A f-k\ , . 

a = -(t0' (29) 

-rxe-*^ = ^{±), (30) 



e x+(t> ( k 

From equations (27) and (31), 



x = • (31) 



A = re^SvrTo - (32) 
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In the light of equation (32), we may make some considerations on the phases of 
graviatational collapse: 

1) In a first phase, A increases very fastly, due to the important flux of infalling 
matter and the insignificancy of pre-Hawking radiation. 

2) When A reaches a certain value, we arrive to a phase of "stability", where the 
infalling flux of collapsing matter is exactly compensated by the outgoing flux of pre- 
Hawking radiation. 

3) Finally, when the "outer layers" of infalling matter have already got exhausted, 
the pre-Hawking term prevails and A will diminish. 

In this paper, we are going to focus our study mainly on the second phase. 



2.2 The stability phase 

In the stability phase, e~ A << 1, so that 

8vrT w (33) 

(the approximations that we are going to perform in this section are not good for small 
values of r). 

The condition for the stability phase consits on imposing X st = 0: 

= re-*- (1) - (34) 
<f>st = -\ st + ln(r 2 ) -In (2k). (35) 



In this situation, 



Consequently, 



Xst = (36) 



Tl = 0, (37) 
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From equation (35), 



r 

From equations (41) and (42), 



'st = -Kt + l*-Kt- (42) 



p^st 

Kt « — , (43) 



(44) 



By integrating, 



r 

d\ st ^ e A ^ 
dr r 

e-^dXst « --. (45) 
r 



- e" Ast » - ln(r), (46) 
A st «-/n(-^i + /n(r)). (47) 



2.3 On contour conditions 



The function (f) has a null value outside the edge of the collapsing star. If we call M the 
total mass of the star, its radius R is given by an expression of the following type [Pin]: 



i? = 2M + e £ w i , (48) 



where M = M(t). 
Thus, 



MR) = -\ s t(R) + In (i? 2 ) - In (2k) = 0, (49) 

X st (R) = In (i? 2 ) - In (2k) » -In {-K x + ln(R)) . (50) 

2k 

Kt nlnR-rjp (51) 
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